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X) ■ Abstract 

OO I The Vincent-Phatak procedure for solving the momentum-space Schrodinger 

equation with combined Coulomb-plus-short-range potentials is extended to 

> . 

CN . angular momentum states coupled by an optical potential — as occurs in spin 

cn 
o 

(N 
O . 

, parameterization is derived and applied to 500 MeV polarized-proton scatter- 



1/2 X 1/2 scattering. A generalization of the Blatt-Biedenharn phase shift 



ing from '^He and ^^C. The requisite high-precision partial- wave expansions 
and integrations are described. 

I. INTRODUCTION 

The theory and equations of quantum mechanics are represented equally well in coordi- 
nate and momentum spaces. Bound states problems, which by definition deal with normal- 
izable wavef unctions, can actually be solved equally well in either space, while scattering 
problems, which in the time-independent Schrodinger theory deal with non-normalizable 
states, are more of a challenge in momentum space. This challenge arises, in part, be- 
cause boundary conditions are more naturally imposed in coordinate space, and, in part, 
because non-normalizable states cannot be Fourier transformed O. In spite of the difiicul- 
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ties, momentum-space calculations are important because fewer approximations are needed 
there to handle the nonlocal potentials arising in many-body and field theories. 

The Coulomb problem in momentum space has actually been "solved" a number of 
times — possibly starting with Fock's study of the hydrogen atom |0 — yet no one numerical 
approach appears to provides the requisite precision for all applications. The real "problem" 
is that the Coulomb potential between a point projectile (P) and a target (T), 

K(k',k) = ^^p(g) , g=|k'-k|, (1) 

has a singularity which must be regularized before a numerical solution is implemented 
0. (The form factor p{q) in (|l|) accounts for the finite size of the target's charge distribution 
and makes the potential well-behaved at large q — but does not remove the singularity at 
g = 0.) Kwon and Tabakin solved the bound-state problem with the potential (|I]) by 
using Lande's technique of subtracting a term from ([^), which makes its integral finite, 
and then adding in a correction integral. Alternatively, Cieply et al. solved the bound- 
state problem by using the Vincent-Phatak (VP) procedure [|^, which deals with the Fourier 
transform of a Coulomb potential which has been cut off beyond some radius Rcut, 

Vr'i^', k) = [p{q) - cos(gi?,„,)] ■ (2) 

While V^'^"*(k', k) is clearly finite as g ^ 0, it produces wavefunctions which must have their 
asymptotic behavior corrected. 

The VP procedure was originally formulated for intermediate-energy pion scattering from 
light nuclei where it provided sufficient accuracy 0. However, the accuracy has become 
a concern for intermediate-energy proton scattering where the proton's much larger mass 
leads to correspondingly larger momentum transfers and correspondingly greater numbers 
of partial waves. Crespo and Tostevin and Picklesimer et al. |T0| have documented 
difficulties with the VP procedure, difficulties which appear as a sensitivity of the computed 
phase shifts to the cutoff radius, or a several-percent error in the phase shift when compared 
to coordinate-space calculations. Both references suggest algorithms to reduce the errors. 
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Alternatively, Elster et al. applied the two-potential formula to the Coulomb and nuclear 
potentials and outlined an approach requiring multiple, numeric Fourier transforms between 
coordinate and momentum spaces. In contrast, Arrellano et al.'s study of intermediate- 
energy proton scattering from spinless nuclei simply made the VP procedure sufficiently 
precise by using some high-precision partial wave expansions developed by Eisenstein and 
Tabakin |jl3|. (As a check, they transformed the potentials to coordinate space and solved 
the equivalent integro-differential equation.) 

In the present paper we generalize the Blatt-Biedenharn phase shift parameterization 
and the VP procedure in order to handle intermediate-energy proton scattering from spin 
1/2 nuclei in which tensor forces couple states of differing angular momenta. In Sec. |I| 



we reformulate the VP procedure for uncoupled channels, in Sec. |TT| we present our new 
formulation for coupled channels, and in Sec. |^ we give some details of an application to 
500 MeV proton scattering from ^He. 



II. UNCOUPLED STATES (0 x 0, x 1/2) 

Vincent and Phatak formulated their procedure in terms of phase shifts. We reformula- 
tion it in terms of T-matrix elements to avoid unnecessary conversions between phase shifts 
and amplitudes. The summed nuclear plus cut-off Coulomb potential V, can be used in the 
Lippmann-Schwinger equation: 

Ti±{k , k) = Vi±{k ,k) + - p dp — — , (3) 

IT Jo E + ie — E[p) 

which is then solved in the normal way since Vi±{k\ k) is well-behaved in momentum space 
(short-ranged in coordinate space). If the r-space nuclear potential vanishes beyond a range 
i?, and if the Coulomb cutoff-radius Rcut is chosen larger than i?, then the wavefunction in 
the intermediate region, R < r < Rcut, can be expressed in terms of a free wave shifted by 
an intermediate phase shift 6i [0: 

Uj=i±i/2{R < r < Rent) = Ne'''^ [sin5;± Gi{kr) + cos5;± Fiikr)] , (4) 
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= N[Fi{kr)+fi±H^+\kr)], (5) 
Ti± = e'''± sm 5i± = -PeTi± , Pe = 2ko \\ \/ , 6 

where 5j=i± is the intermediate phase shift arising from the short-range potentials and is 
a normahzation constant. Here I is the orbital angular momentum, j = / ± 1/2 is the total 
angular momentum, and we have used two equivalents forms for the wavefunctions [Q. 

The wavefunction for r > Rcut can be expressed in terms of a phase shift 5i (the amount a 
point-Coulomb wave is shifted), which, in turn, is determined by matching the intermediate 
wavefunction to an outer one at r = Rcut- The outer wavefunction for r > Rcut has the 
same form as the intermediate one, but with the free waves replaced by Coulomb waves and 
the intermediate phase shift 5i± replaced by the phase shift relative to point-Coulomb waves 

ui± (r > R,^t) = N' [Fi {v,kr) + H^^^ (r/, kr)] , f 4 = e*^'± sin . (7) 

Here TJj_ is the T matrix for scattering from short range forces in the presence of a point- 
Coulomb force, 1] = ZpZtC^/v is the Sommerfeld parameter, Fi{ri, kr) is the regular Coulomb 
function, and Hl~^\r], kr) is an the outgoing wavefunction; for example, 

Hj^{r], kr) = Gi^i], kr) ± iFi(r], kr) ~ exp{±i[kr — ln/2 + ai — r]\n{2kr)]} . (8) 

The requirement that the logarithmic derivative of ui±{r < Rcut) match ui±{r > Rcut) at 
r = Rcut determines T^'^ as: 



[F,,H^^\v)]+fiAH;^\Hl^\v)] 
where the brackets indicate Wronskians 



, m r7T(+) rri+)f - ' 



Finally, the non-spin-flip scattering amplitude can be expressed in terms of the outer, 



Coulomb-modified phase shifts 6f ||I4 



fid) = UO) + ne) , (10) 

UO) = -^^-^(^exp{2z [o-o - r/lnsin(e/2)]} , (11) 

1 oo 

r'(^) = ^ E(2^ + (e'^'' - l) ^Kcos^) . (12) 



Here /^^ is the amplitude for scattering from a point-Coulomb potential and /"''^ is the am- 
plitude for scattering from the short-ranged potentials in the presence of the point Coulomb 
force. 



III. COUPLED STATES (i x i) 

If the interaction couples angular-momentum states, we must generalize the VP 
method — even if the Coulomb force does not directly couple the states. When two non- 
identical spin 1/2 particles interact through a tensor force, the total angular momentum j 
remains a good quantum number yet there is coupling within the triplet spin state as well as 
between the triplet and singlet states. Accordingly, is generalized to the coupled integral 
equations: 

where the sum is over the coupled states. If the two particles interact through an optical 
potential, the phase shifts are complex and the S matrix non-symmetric. In the Appendix 
we generalize the conventional Blatt-Biedenharn NN parameterization ||l5| , p!6| , |r7| and show 
that the S matrix elements have the form: 

Si'=j±i,i=j±iiko) = S±± = 5w - 2ipETii>{ko, ko) (14) 

S^^ = (cos e+_ cos e_+e^'^++ + sin e+_ sin e.+e^*"^— ) / detU , (15) 

= (sin e_+ cos e.+e^'^++ - cos e_+ sin e.+e^*^" )/ det U , (16) 

= (sin e+_ cos e+_e^*'^++ - cos e+_ sin e+.e^*^" )/ det U , (17) 

S__ = (cose+_cose_+e^^''— +sine+„sine_+e2^'^++)/detf/ , (18) 

det U = cos e+_ cos e_+ + sin e+_ sin e_+ . (19) 

To apply the VP procedure to channels coupled by an optical potential, we 1) transform 
to a new basis in which there is no channel coupling, 2) match the interior wavefunction 
to point-Coulomb ones in this basis for which the S-matrix is diagonal, and then 3) return 



to the original basis to calculate the scattering observables. One implementation of these 
steps would be to take our S matrix elements computed via (|T3|) and ([T^), assume they have 
the forms (|T5|)-(|1^) in terms of phase shifts and coupling parameters, and then search for 

the {6 , 5++, e+_, e_+) which satisfy these transcendental equations. The 5's would be the 

phase shifts in the basis in which S is diagonal — even though we never explicitly transform 
to that basis. 

The implementation we have used is more direct and self-testing. We compute the 
solution of ([13|) , form the nondiagonal S matrix, 

S — u S — 



[S] 



(20) 



and then explicitly diagonalize it with the similarity transformation: 

[s'] = [u][s][ur " ° 







[U] 



1 

S- + 



S+- 



A+-5++ 



\--S- 
1 



[U] 



-1 



A_-S- 
1 



detU 



A 



± 



SV + AS+^S, 



(21) 

(22) 
(23) 



We now effectively deduce the intermediate 6 phase shift from the diagonal elements, do 
the VP matching with the corresponding T^s as in (|^) (this effectively determines the final 
phase shift 6'^), and then we use the original U matrix to transform back to the basis in 
which we calculate observables: 



[5-] = [U]-^[S'][U] . 



(24) 



Even though the method is guaranteed to diagonalize the S matrix, as an internal test we 



check that \Su?\ < 1 and that 



Sill 



< 1. 



IV. COMPUTATIONAL DETAILS 



We have modified the LPOTp code |T3] to include the Coulomb potential in the different 
spin channels. As a first test of our precision we solved the Lippmann-Schwinger equation 
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with a point Coulomb potential and checked that our answers reproduced the point Coulomb 
phase shifts ai [after removal of the riln{2kRcut) term in (^)]. We concluded from this severe 
test that 48-64 grid points are required to solve the Lippmann-Schwinger equation and 
obtain four-five place precision in ai [there is enough cancellation that three-place precision 
does not reproduce the point Coulomb scattering amplitude (O)]- After an overall e^*°"° 
is factored out from the sum in (|Tll), good agreement for was found (indistinguishable 



from the analytic amplitude on a five-decade semi- logarithmic plot). 

As the next test we computed pure-Coulomb scattering of 415 MeV protons from the 
charge distribution of ^He. We were able to obtain essentially perfect reproduction of the 
Born-approximation amplitude, 

f finite ^ fot^Q^p^q^ ^ (25) 

which proves that we can include short-ranged effects — in addition to the long-range 
Coulomb force — with precision of at least 0{a^). To actually obtain this agreement we 
used 48 grid points in the solution of the Lippmann-Schwinger equation (|TB|), and increased 
the precision of our partial-wave projection: 

Vi{k\k) = j\ik',k)Piicose',,)d{cosei,) , (26) 

until the partial wave summation, 

1 ^max 

V{k',k) ^ — Y^{2l+l)Vi{k',k)Pi{cosek>k) , (27) 

reproduced all oscillations present in V{k', k). We show a reproduction of this type in Fig- 
ure where the many oscillations arising from the cos{qRcut) term in the cut-off Coulomb 
potential (H) is evident. We obtained six-place reproduction of V(k', k) using l^ax = 48 
partial waves and 96 integration points in the partial- wave projection (^). Ten-place repro- 
duction demanded Imax = 96. We expect these number to scale as kR, and so larger nuclei 
or higher energies will require more partial waves and grid points. For these calculations 
we used analytic nuclear form factors [l^, though we also were successful for ^^C using 



numerical Fourier transforms of Wood-Saxon densities [|T9[. However, noise and instability 
do appear for form factors which fall off slowly in q. 

An important requirement on the VP method is that the matching radius, which we take 
equal to Rcut-, be larger than the range of the nuclear force (in order to be able to express the 
outer wavefunction as a linear combination of Coulomb waves). However, increasing R^ut 
makes the cut-off Coulomb potential more oscillatory and more difficult to reproduce. In fact, 
it was the sensitivity to changes in Rcut which led Ref. to search for an alternative to the 
matching method. We find that using Rcut < 5 fm produces unstable results (presumably 
cutting off the nuclear potential), but, as seen in Figure |^, we obtain stable results for 
6 fm < Rcut < 10 fm. 

In Figure ^ we compare the nuclear-plus Coulomb cross section and polarization (solid 
curves) to those calculated without Coulomb (dashed curves). The exact handling of the 
Coulomb potential is seen to have a significant, although small, effect in the semilog plot 
of da/dQ, and a more pronounced effect for Aqq^q. Not plotted, because they essentially 
overlap the exact results, are ones in which the Coulomb potential is handled in impulse 
approximation: 

m ^ ud)p{q) + ne) , (28) 

with f^{d) the scattering amplitude for pure nuclear scattering. 



V. CONCLUSION 

We have extended the Vincent- Phatak procedure for the exact inclusion of the Coulomb 
potential in momentum space to calculations of proton scattering from spin 1/2 nuclei in 
which spin-dependent forces couple angular-momenta states. As part of that extension we 
also generalized the Blatt-Biedenharn phase shift analysis for the scattering of two spin 1/2 
particles to cases where the S matrix is no longer symmetric. Although our formulation 
and calculational procedure is for a more complicated spin case, we confirm the finding of 
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Arrellano et al. [|T2[ that the VP procedure can be made sufficiently accurate for intermediate- 
energy proton scattering if high-precision partial-wave expansions and large numbers of 
partial waves are used. 
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APPENDIX: GENERALIZATION OF BLATT-BIEDENHARN CONVENTION 



The conventional phase shift analysis must be extended when the angular momentum 
channels are coupled. Blatt and Biedenharn did this first for the mixing of the / = j ± 1 
states within the nucleon-nucleon triplet state ||l6| , p!7| by assuming that the mixed states 
have the asymptotic forms: 

lim u, i=i-i(r) = A , e-'[^''-(^-^)"/2] _ ^ ilkr-{j-i)n/2] .^^^ 

r— >oo J' J ^ ' ' 

The S matrix is defined by the relation among the A's and B's: 



(A3) 



B+ _ S++ 5+- 

For NN scattering below pion production threshold, S must be unitary because flux is 
conserved, and symmetric because all terms in the Schrodinger equation are real. For that 
case, the most general form for S, a unitary and symmetric 2x2 matrix, is given by a 
similarity transformation with mixing parameter e, 



[S] = [U]-\e'^%U] 



[U] 



cos ej sm 



sm ej cos tj 



2i(54 



(A4) 
(A5) 
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where 5++ = 5i=j+i^i'=j+i and S — = Si=j-i^i'=j^i. 

When deahng with an optical potential, the S matrix is no longer unitary — which means 
the phases shifts become complex, as well as no longer symmetric — which means there are 
now two mixing parameters. We assume ( [A^ ) to be valid with the more general transfor- 
mation matrix: 



[U] 



cosej 



sm 6- 



, [U] 



-1 



1 



-sme+_ cose i 

det U = cos e+_ cos e_+ + sin e+_ sin e_ 



detU 



cos e ^ — sm 

sin e+_ cos e+_ 



(A6) 
(A7) 



This leads to the S matrix elements given in ([T5| ) - ([l8|) which reduce to the standard, coupled 
p!5| , p!6| if e+_ = e_+, and to the standard uncoupled case if e+_ = e_+ = 0. Stapp ||T5[ 



case 



also gave a parameterization of the S matrix in terms of the "bar" phase shifts which are, 
in some cases, more convenient in the parameterization of data. These bar phases, however, 
are not the ones introduced here, and, in fact, do not provide a diagonal representation of 
the S matrix. 
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FIGURES 

FIG. 1. The nuclear plus Coulomb potentials in momentum space for the spin triplet state with 
nis = rUg' = 1 as a function of the cosine of the angle between k and k'. The summation ([27| ) of 
partial-wave potentials essentially overlaps the input potential. 

FIG. 2. The differential cross section for 500 MeV proton scattering from '^He. Calculations 
performed using a cutoff radius in the range 6fm < Rcut ^ 10 fm fall within the two curves. The 
experimental data are from Hausser et al. pO| . 

FIG. 3. The differential cross section and analyzing power (unpolarized target, projectile po- 
larized in normal direction) for 500 MeV proton scattering from '^He. The solid curves gives the 
exact results using the VP method and the dashed curves gives the results if no Coulomb force is 
included. The experimental data are from Hausser et al. pO|l . 
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